An  eddy-resolving  ocean  reanalysis  using 
the  1/12°  global  HYbrid  Coordinate  Ocean 
Model  (HYCOM)  and  the  Navy  Coupled 
Ocean  Data  Assimilation  (NCODA)  scheme 


O.M  Smedstad1,  EJ.  Metzger2,  R.A.  Allard2,  R.  Broome1, 

D.S.  Franklin1  and  AJ.  Wallcraft2 

1QinetiQ  North  America 
2Naval  Research  Laboratory 


Layered  Ocean  Model  Workshop 
21-23  May  2013 
Ann  Arbor,  Michigan 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

MAY  2013  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2013  to  00-00-2013 

4.  TITLE  AND  SUBTITLE 

An  eddy-resolving  ocean  reanalysis  using  the  l/12degree  global  HYbrid 
Coordinate  Ocean  Model  (HYCOM)  and  the  Navy  Coupled  Ocean  Data 
Assimilation  (NCODA)  scheme 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

Naval  Research  Laboratory, 1  Stennis  Space  Ctr  ,Stennis  Space 

Center, MS, 39529 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS  (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ARSTRATT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  3S 

unclassified  unclassified  unclassified  Report  (SAR) 

24 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


HYCOM/NCODA  Ocean  Reanalysis 

Of  those  ocean  reanalyses  performed  to  date,  only  a  few 
have  eddy-permitting  resolution  and  none  are  capable  of 
fully  resolving  oceanic  mesoscale  features  (eddies  and 
current  meanders)  across  the  globe 

This  project  addresses  the  need  for  a  long  time  period 
eddy-resolving  ocean  reanalysis 

Funded  by  the  DoD  Modeling  and  Simulation 
Coordination  Office  (M&S  CO) 

Goal  to  the  sponsor:  provide  physically  consistent 
environmental  scenarios  for  planning,  wargaming  and 
scenarios  to  support  the  warfighter 

Numerous  other  applications  and  research  opportunities 
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HYbrid  Coordinate  Ocean  Model 

Tri-pole  latitudinal  grid  resolution  (km) 


Curvi-linear 
grid:  north 
of 47°N 


Mercator 
projection: 
66°S  to 
47°N 

Uniform 
cylindrical: 
south  of 
66°S 


•  32  hybrid  coordinate  surfaces,  thermobaricity,  o2* 

•  K-Profile  Parameterization  (KPP)  mixed  layer  model 

•  Monthly  river  runoff 

•  Surface  salinity  relaxation  to  U.S.  Navy  GDEM4  climatology 

•  Thermodynamic  "energy  loan"  ice  model 


Atmospheric  Forcing 


NCEP  Climate  Forecast  System  Reanalysis  (CFSR) 

Surface  Specific  Humidity  (kg/kg  x  102) 


•  Time  frame:  1993-2012  (altimeter  period) 

•  Horizontal  resolution:  0.3125°  gaussian 
•Temporal  resolution:  1-hourly 

•  QuikSCAT  scaling 

•  Inputs: 

•  Bulk-derived  wind  stress 

•  Wind  speed 

•  Radiative  fluxes 

•  Thermal  fluxes 

•  Precipitation 


Net  Surface  Shortwave  (W/m2) 


CFSR-sea  so  rad  20040321  ooz 


Precipitation  (m/s  x  106 


ChSK-std  Precip  2004032T00Z  m/s  +  1.0 L6 


Spin-up  and  Reanalysis  Simulations 

•  Spin-up  1/12°  non-assimilative  global  HYCOM 
with  CFSR  climatology  (12  model  years) 

•  Extend  climatological  spin-up  with  1993-2012 
1-hourly  CFSR  forcing  using  1/12°  non- 
assimilative  global  HYCOM 

•  Extend  climatological  spin-up  with  Oct  1992- 
2012  1-hourly  CFSR  forcing  using  1/12° 

assimilative  HYCOM/NCODA 
—  Begin  in  Oct  1992,  currently  in  May  2000 


Navy  Coupled  Ocean  Data  Assimilation 

(NCODA) 


Ocean/Ice 

Obs. 


SST:  G AC/LAC  MCSST, 

GOES,  Ship,  Buoy 

Profiles:  XBTs,  CTDs,  Argo 
floats 

Buoys:  Fixed  and  drifting 
Altimeter  SSH  anomaly 
SSM/I  sea  ice  concentration 


Sequential  Incremental  Analysis  Update  (IAU) 
Analysis-Forecast-Analysis 
NCODA  analysis  once  per  day  that  is 
incrementally  inserted  over  a  6  hour  window 


Ocean  QC 


Innovations 


3DVar 


Increments 


3DVar  -  simultaneous  analysis  of 
ice  concentration  and  5  ocean 
variables:  temperature,  salinity, 
geopotential,  layer  pressure,  and 
velocity  (u,v) 


First  Guess 


HYCOM 
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Altimeter  Observations 


January  1, 1993 


August  26, 1999 


February  19, 2000 
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Observations  to  be  Assimilated  via  NCODA 

August  16  -  September  16  1999 


SSH  Observations:  3-day  data  window  Modular  Ocean  Data  Assimilation  System  (MODAS) 
SST  Observations:  24  hour  data  window  used  as  vertical  projection  of  the  satellite 
Profiles:  12-day  data  window  altimeter  observations 
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Mean  Bias  (C)  RMS  Error  (C) 


Temperature  Verification  in  the  Pacific 

August  1998  to  August  1999 
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SST  Mean  Error 


1999 


SSH  Variability 


191  SSH  Variability  -  Year  1993 1999 


0°  _ 90°E _ 180! _ 90°W _  0° 


1993-1999 


Oct  92  -  May  07  SSH 
variability  based  on 
T/P,  ERS-1  and  ERS-2 
altimeters  (from 
Collecte, 
Localisation, 
Satellites  (CLS)) 
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Mean  SSH 

1993-1999 


Gulf  Stream  region  Kuroshio  region 
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Kuroshio  and  Gulf  Stream  SSH  with  SST- 

based  frontal  analysis 


SSH  May  16  2000 


SSH  April  25  2000 


sea  surf,  height  20000516  [19.1] 
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Frontal  analysis  <  4  days  old  =  white. 

Frontal  analysis  >  4  days  old  =  black 

Frontal  analysis  performed  by  the  Naval  Oceanographic  Office 


Mean  SSH 

Gulf  of  Mexico  Region 


Gulf  of  Mexico  SSH  and  SST  with  SST- 

based  frontal  analysis 

March  3  1999 


Frontal  analysis  <  4  days  old  =  white, 
analysis  >  4  days  old  =  black 

Frontal  analysis  performed  by  the  Naval  Oceanographic  Office 
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1993-1999  Mean  SSH  in  Luzon  Strait 


Forced  with  1  hourly  CFSR  winds  and  thermal  forcing 


Velocity  Cross-section  Across  Luzon  Strait 

Sb-ADCP  data  (top)  vs.  1/12°  Global  HYCOM  (bottom)  in  the  upper  300  m 
Section  along  120.75°E  between  Taiwan  and  Luzon 

19°N  20°N  21°N  22°N 


Westward  velocity: 
40+  cm/s 


Westward  velocity: 
30+  cm/s 


Sb-ADCP  data  from  Liang  et  al.  (2003, DSR  Pt.  II),  1991-2000 
1993-1999  mean  from  HYCOM  forced  with  1  hourly  CFSR  winds  and  thermal  forcing 


Velocity  Cross-section  Along  Luzon  Strait 


Sb-ADCP  data  (top)  vs.  1/12°  Global  HYCOM  (bottom)  in  the  upper  300  m 
Section  along  21°N  between  118. 5°E  and  124.0°E 

119°E  120°E  121°E  122°E  123°E  124°E 


Northward  velocity 
of  western  core: 
60+  cm/s 


Northward  velocity 
of  western  core: 
40+  cm/s 


Sb-ADCP  data  from  Liang  et  al.  (2003,  DSR  Pt.  II),  1991-2000 
1993-1999  mean  from  HYCOM  forced  with  1  hourly  CFSR  winds  and  thermal  forcing 


Velocity  Cross-sections  East  of  Taiwan 

Sb-ADCP  data  (top)  vs.  1/12°  Pacific  HYCOM  (bottom)  in  the  upper  300  m 

Sections  at  22°N,  23°N,  24°N  and  25°N 

22°N  23°N  24°N  25°N 


121°E  122°E  123°E  121°E  122°E  123°E  121°E  122°E  123°E  121°E  122°E  123°E 


Sb-ADCP  data  from  Liang  et  al.  (2003,  DSR  Pt.  II) 

1993-1999  mean  from  HYCOM  forced  with  1  hourly  CFSR  winds  and  thermal  forcing 


Note  how  the  two-core  Kuroshio  merges  to  a  single  jet  in  both  the  observations 
and  HYCOM  from  the  south  to  north  along  the  Taiwan  coast 


WOCE  PCM-1  Transport 


From  Johns  et  al.  2001 
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WOCE  PCM-1  Transport 


WOCE  PCM-1  transport  (black)  via 
the  adjusted  geostrophic  method 
from  Johns  et  al.  (2001)  versus  full 
water  column  transport  from  the 
global  HYCOM/NCODA  Ocean 
Reanalysis  (red)  using  those  model 
gridpoints  between  moorings  Ml  - 
M4.  A  ten  day  filter  has  been 
applied  to  both  time  series.  The 
unfiltered  mean  and  standard 
deviation  for  PCM-1  are  22.0  ±  4.8 
Sv  while  that  for  global  HYCOM  is 
23.7  ±  9.8  Sv.  The  correlation 
coefficient  is  0.78. 
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Computational  Requirements 

•  Computer  time  via  the  DoD  High  Performance 
Computing  Modernization  Office 

•  Currently  integrating  the  ocean  reanalysis  on 
the  Navy  DoD  Supercomputing  Resource  Center 
(DSRC)  IBM  iDataPlex 

•  Using  949  processors 

-  Integrate  up  to  22  model  days  every  24  hrs  of  wall 
time 

-  It  will  take  ~8  months  to  integrate  the  remainder 
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Output  and  Storage 

•  HYCOM  3D  native  grid  archive  files  (compressed): 

-  Single  hour:  ~7  Gb 

-  Saving  3-hourly  output: 

•  ~20  Tb  /  model  year 

•  ~340  Tb  for  the  entire  reanalysis 

•  HYCOM  3D  constant  .08°  grid  (±80°  lat)  netCDF  files 
remapped  to  40  z-levels  (compressed): 

-  Single  hour:  ~1.2  Gb 

-  Saving  3-hourly  output: 

•  ~3.5  Tb  /  model  year 

•  ~59  Tb  for  the  entire  reanalysis 

•  Subset  of  the  output  will  be  placed  on  the  hycom.org 
data  server 

-  When  is  still  to  be  determined 


24 


